Casimir interaction between a plate and a cylinder 
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We find the exact Casimir force between a plate and a cylinder, a geometry intermediate between 
parallel plates, where the force is known exactly, and the plate-sphere, where it is known at large 
separations. The force has an unexpectedly weak decay ~ L / (H 3 ln(H / R)) at large plate-cylinder 
separations H (L and R are the cylinder length and radius), due to transverse magnetic modes. 
Path integral quantization with a partial wave expansion additionally gives a qualitative difference 
for the density of states of electric and magnetic modes, and corrections at finite temperatures. 
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With recent advances in the fabrication of electronic 
and mechanical systems on the nanometer scale quan- 
tum effects like Casimir forces have become increasingly 
important 0,0- These systems can probe mechanical os- 
cillation modes of quasi one-dimensional structures such 
as nano wires or carbon nanotubes with high precision Q • 
However, thorough theoretical investigations of Casimir 
forces are to date limited to "closed" geometries such as 
parallel plates Q or, recently, a rectilinear "piston" 
where the zero point fluctuations are not diffracted into 
regions which are inaccessible to classical rays. A notable 
exception is the original work by Casimir and Polder on 
the interaction between a plate and an atom (sphere) at 
asymptotically large separation 

In this Letter we consider the electrodynamic Casimir 
interaction between a plate and a parallel cylinder (or 
"wire"), both assumed to be perfect metals (see inset of 
Fig. 2J. We show that the Casimir interaction can be 
computed without approximation for this geometry. We 
believe that the methods presented here may yield exact 
solutions for other interesting geometries as well. This 
geometry is also of recent experimental interest: Keeping 
two plates parallel has proved very difficult. The sphere 
and plate configuration avoids this problem, but the force 
is not extensive. The cylinder is easier to hold parallel 
and the force is extensive in its length 0j. 

Casimir interactions, while attractive for perfect met- 
als in all known cases, depend strongly on geometry 
Consider the Casimir interaction energy (discarding sep- 
aration independent terms) at asymptotically large H 
for three fundamental geometries which differ in the co- 
dimension of the surfaces 0: two plates, plate-cylinder, 
and finally, plate-sphere, corresponding to co-dimension 
1, 2, and 3, respectively. It is instructive to consider 
both a scalar field which vanishes on the surfaces (D = 
Dirichlet) and the electromagnetic field (EM). For par- 
allel plates (area A) E ~ hcA/H 3 in both cases 0. For 



a plate and a sphere of radius R, E ~ hcR/H 2 for 
the Dirichlet case, as compared with E ~ hcR 3 /H 4 for 
the EM case @. Based on these results, expectations 
for the plate and cylinder geometry might range from 
~ hcLR 2 j H 4 , proportional to the cylinder volume, to 
~ HcLR/H 3 , proportional to its surface area, or even 
~ HcL/H 2 with a potential non-power law dependence 
on the radius. 

A simple but uncontrolled method for study of non- 
planar geometries is the proximity force approximation 
(PEA) , where the system is treated as a sum of irffinites- 
mal parallel plates Applied to the plate-cylinder 

geometry, the PFA yields .Bp fa = — ^qTt 3 hcL ^ R/ 2a 5 
to leading order in a/R, where a = H — R. Other 
approximations include semi-classical methods based on 
the Gutzwiller trace formula |lll| . and a recent optical 
approach which sums also over closed but non-periodic 
paths 01 ■ F° r large separations, a multiple scattering 
approach is available |l3j . but has not been adapted to 
this geometry. For the Dirichlet case, a Monte Carlo ap- 
proach based on worldline techniques has been applied 
to the plate-cylinder case [lij . 

Our result provides a test for the validity of these ap- 
proximate schemes, and also provides insight into the 
large distance limit. In particular, we find the unex- 
pected result that the electrodynamic Casimir force for 
the plate-cylinder geometry has the weakest of the pos- 
sible decays, 



F = 



hcL 



8tt H 3 \n(H/R) 



(1) 



as H/R — > oo. The same asymptotic result applies to a 
scalar field with Dirichlet boundary conditions. Interest- 
ingly, the decay exponent of the force is not monotonic 
in the number of co-dimensions: (—4, —3 — e, —5) for co- 
dimension (1,2,3) respectively. In contrast the Dirichlet 
case is monotonic with exponents (— 4,— 3 — e, —3). 
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In the remainder we derive these results, summarized 
in Eqs. JSJ-lJHJl, using path integral techniques. Our ap- 
proach also yields the distance dependent part of the 
density of states, which contains the complete geometry 
dependent information of the photon spectrum, and is 
useful for computing thermal contributions to the force. 

The translational symmetry along the cylinder axis en- 
ables a decomposition of the EM field into transverse 
magnetic (TM) and electric (TE) modes 0] which arc 
described by a scalar field obeying Dirichlet (D) or Neu- 
mann (N) boundary conditions respectively. We can 
compare our TM results to recent Monte Carlo Dirich- 
let results 0| . Moreover, the mode decomposition turns 
out to be useful also in identifying the physical mecha- 
nism behind the weak decay of the force, which at large 
distance is fully dominated by D modes. 

Our starting point is a path integral representation 0] 
for the effective action whic h y ields a trace formula for 
the density of states (DOS) |l7]|. The latter is then eval- 
uated using a partial wave expansion. The DOS on the 
imaginary frequency axis is related to a Green's function 
by p(iqo) = (2<7o/tt) / d 3 xG(x, x; q ), where G(x,x';q ) 
is the Green's function for the scalar field with action 
S=±J d 3 x(\V^\ 2 + (q°) 2 $ 2 ). The effect of boundaries 
on the Green function can be obtained by placing func- 
tional delta functions on the boundary surfaces in the 
functional integral^]. By integrating out both the field 
$ and the auxiliary fields which represent the delta func- 
tions on the surfaces, one obtains the trace formula 0| 



can be expressed as, 



E 



hcL 
8^2 



dqodqi In 



dct A/(g ,gi) 
det Moo(qo,qi) 



(3) 



The elements of the matrix M are labeled by the in- 
teger index m = — oo, . . . , oo coming from the compact 
^-dimension of the cylinder, and the momentum q 2 along 
the other direction parallel to the plate, to read 



M 



^[m,m'] ^[m,q' 2 ] 
B [q 2 ,ra'\ ^[l2,q' 2 ] 



(4) 



The matrix A 
A, , = A 

m,m — -^n 



m,m'l is diagonal, with elements 
i = I m (ru)K m (ru) for D and 
Am = (u/H) 2 I m (ru)K' m (ru) for N modes. The ma- 
trix C also has only diagonal elements C[q 2 ,q 2 ] = 

C(q 2 ) = H/(2y / u 2 + u 2 ) for D and C(q 2 ) = 

— y / u 2 + u 2 / (2H) for N modes. The off-diagonal ma- 



trix B is non-diagonal with B\ 



m,Q2\ 



B m (q2) = 



TtHe-^ u2+u *I m {ru)[u/{u 2 + V u 2 + u 2 2 )] m / V u 2 + u\ for 

D and B m (g 2 ) = (nu/ ' H)e-^ u2 ^I' m (ru)[u/ '(u 2 + 
y^u 2 + u 2 )]~ m for N modes. Here, we have defined the 
dimensionless combinations u = H y q$ + qf, u 2 = Hq 2 
and r = R/H. The determinant can be obtained 
straightforwardly, and the total energy can be decom- 
posed to the sum of D and N mode contributions, as 



N , 



(5) 



with 



M9o) = --^-Trln(A/A/ 00 1 ) 
7T dq 



(2) 



where Sp(qo) is the change in the DOS caused by mov- 
ing the plate and cylinder in from infinity. The infor- 
mation about geometry is contained in the matrix M 
of the quadratic action for the auxiliary fields, given by 
M aj3 (u, u';q ) = G (s Q (u) - sp(u');q ) for D and by 
M a p(u,u';q ) = d nm ^d np ^G (s a (u) - sp(u');q ) for 
N boundary conditions; Go = e~ 9o ' x '/47r|x| is the free 
space Green function, d na is its derivative normal to the 
surface, and s a (u) parametrizes the surfaces (which are 
numbered by a = 1, 2) in terms of surface coordinates 
u. M^ 1 is the functional inverse of M at infinite sur- 
face separation. The trace in Eq. J5J) runs over u and 
a. For the cylinder with its axis oriented along the x\ 
direction we set Si(xi,(p) = (x\, Rsm((p), Rcos(ip)) and 
for the plate s 2 {xi 1 x 2 ) = (xi,x 2 ,H) (see inset of Fig. [Q. 

The Casimir energy of interaction is given by E = 
(he/ 2) J Q dqoqo5p(qo). After transforming to momen- 
tum space, M, the Fourier transform of the matrix M 
has block diagonal form with respect to qo and the mo- 
mentum qi along the cylinder axis, so the Casimir energy 



$ A (r) 



4tt Jo 



duu\n [dct(l - N x (u,r))] . (6) 



The matrix N x (u) is given in terms of Bessel functions, 



K„(u,r) 



I v (ru) 
K^(ru) 



K fi+V (2u), 



for D modes and 



N N 



(u, r) 



K(ru) 
K'(ru) 



K„ +V (2u) 



(7) 



(8) 



for N modes. The determinant in Eq. © is taken with 
respect to the integer indices /Lt, v = — oo, . . . , oo. If the 
matrix N x is restricted to dimension (21 + 1) with TV^j 
as the central element, it then describes the contribution 
from I partial waves, beginning with s-waves for / = 0. 

From Eq. ©, one can easily extract the asymptotic 
large distance behavior of the energy for r = R/H <C 1. 
For Dirichlet modes s-waves dominate, while for Neu- 
mann modes both s- and p- waves (1 = 1) contribute at 
leading order in r. The two cases differ qualitatively, with 



$>)=--L-L, and d>» 
167T In r 



32tt 



(9) 
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For H ^> R the EM Casimir interaction is dominated 
by the D (TM) modes. Note that a naive application of 
the PFA for small r, where it is not justified, yields the 
incorrect scalings $ £> (r) = $ JV (r) ~ r. 

The natural expectation from the Casimir-Polder re- 
sult for the plate-sphere interaction, that the force at 
large distance is proportional to the volume of the cylin- 
der, is incorrect. The physical reason for this difference 
is explained by considering spontaneous charge fluctua- 
tions. On a sphere, the positive and negative charges 
can be separated by at most distances of order R <C 
H . The retarded van der Waals interactions between 
these dipoles and their images on the plate leads to the 
Casimir-Polder interaction [6|. In the cylinder, fluctua- 
tions of charge along the axis of the cylinder can create 
arbitrary large positively (or negatively) charged regions. 
The retarded interaction of these charges (not dipoles) 
with their images gives the dominant term of the Casimir 
force. This interpretation is consistent with the difference 
between the two types of modes, since for N modes such 
charge modulations cannot occur due to the absence of 
an electric field along the cylinder axis. Eventually, for a 
finite cylinder, in the very far region H 3> L, the charge 
fluctuations can be considered again as small dipoles, and 
the Casimir-Polder law is expected to reappear, making 
the force proportional to the volume of the cylinder LR 2 . 

We next consider arbitrary separations, and use Eq. © 
to obtain the contribution from higher order partial 
waves. A numerical evaluation of the determinant is 
straightforward, and we find that down to even small 
separations of a/R =0.1 the energy converges at order 
Z = 25, whereas for a/R > 1 convergence is achieved for 
I = 4. Fig.njshows our results for Dirichlet and Neumann 
modes and for their sum which is the EM Casimir energy, 
all scaled by the corresponding £"pfa given above |18| . 
Both types of modes show a strong deviation from the 
PFA for a/R > 1, especially the Dirichlet energy. Fig.Q] 
shows also very recent wordlinc-based Monte Carlo re- 
sults for the Dirichlet case at moderate separations [l9j . 
which agree nicely with our exact results. 

Eq. (0 indicates that the Dirichlet dominated force 
vanishes logarithmically as R — > at fixed H . A similar 
result is obtained when the cylinder is replaced by an in- 
finitesimal thin wire, but an UV cutoff is introduced to 
control short wavelength modes. Both results are a con- 
sequence of the fact that the asymptotic form of Eq. @ 
is independent of the actual shape of the cross section 
of the wire, and the cutoff R can be identified with any 
typical scale of the cross section. The leading asymp- 
totic term in Eq. © is also obtained Jg] from the s-wave 
scattering amplitude for the 2-dimcnsional problem of a 
strongly repulsive potential concentrated on the wire. 

The difference between the D and N modes also ap- 
pears in the density of states, which in turn affects the 
temperature dependence of the Casimir force. From 
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FIG. 1: Ratio of the exact Casimir energy to the PFA for 
the geometry shown in the inset. All curves are obtained at 
order 25 of the partial wave expansion, and the accuracy lies 
within the line thickness even at small a/R. The Dirichlet 
data points are from Ref. E3. 



Eq. J2J) we obtain the expression 



q Q HL 



°du, 
— 7 In 

u 



det 



l-N x (^/u 2 + (q H) 2 ) 



det[l- N x (q H)} 



(10) 

which is convenient both for numerical and analytic com- 
putations. Numerical evaluation yields the results shown 
in Fig. |21 for general values of R. Analytical results in 
the limit of small R/H are obtained by considering only 
the s-waves for Dirichlet modes, and the s- and p-waves 
for Neumann modes. For Dirichlet modes we expand in 
1/ ln^oi?), whereas for Neumann modes the small pa- 
rameter is r = R/H. To leading order we find 

SMQo) = + 6q Q H)e- 2c > oH r 2 + G(r 3 ).(llb) 

Fig. fallows for an assessment of the validity range of the 
expansions of Eq. (|llf> which are shown as solid curves. 

These results for the DOS allow us to evaluate for the 
first time finite temperature contributions to the Casimir 
interaction in an open geometry. The difference between 
the free energ y J- and the Casimir energy at T = can 
be written as [13J (kg is Boltzmann constant) 



5T = T - E = irk R T 



dqog(qo)5p(q ), (12) 



with the function g(qo) = J2'kLi s ' m ('2 7T kqoAT)/(iTk) and 
A T = hc/{2nk B T). In the limit R < (H, At) but for 



general H/Xt, we can use the expansion of Eq. to 
obtain to leading order in l/ln(i?/A-r) and R/H for D 
and N modes, respectively, the thermal contributions 



6Fd 



STn 



L 



8 \n(R/X T )H 
k B T LX T R 2 







coth (£) 









64 H A 
7(H/X T ) 2 
sinh 2 (H/X T ) 



7— coth 
At 



6 



H 



H 

Xt . 
cosh(i7/A T ) 



(13a) 
(13b) 
- 20 



sinh J (i7/A T ) 

It is interesting to note that <5.Fn has a minimum at 
H/ Xt = 2.915 . . ., where the corresponding thermal force 
changes from repulsive at small H to attractive at large 
H. At low temperatures, the finite T contributions to 
the Casimir force SF = —dST/OH, 



SF B = ^-k B T 



SF N 



45 

64tt 5 
945 



k B T 
he 

k B T 

he 



HL 



ln(i?/A T ) 

5 

R 2 HL , 



(14a) 



(14b) 



have to be added to Eq. Q for if <C Xt- At larger 
temperatures with C At < -ff, one has the scalings 

SF B ~ k B TLH-' 2 /\n{R/X T ) and SF N k B TLR 2 H~ A . 

At the extreme high temperature limit of At "C R-, only 
thermal fluctuations remain, and ft should disappear from 
the equations. This 'classical limit' is well known for 
parallel plates [2(| and is obtained for smooth, arbit rary 
geometries within the multiple scattering approach [13J, 
and the optical approximation [2l| . (Note that for the D 
modes a subleading ft still survives in the logarithm.) 

Finally, we note that our approach can be extended 
also to multiple wires and distorted beams. Our re- 
sults should be relevant to nano-systems composed of 
1-dimcnsional structures and also to other types of fields 
as, e.g., thermal order parameter fluctuations. 
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